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PREFACE 


The  research  reported  in  this  document  was  conducted  by  Richard  Collins,  Ph.D.,  Biodynamics 
International,  under  a  subcontract  to  ManTech  Environmental  Technology,  Inc.,  in  support  of  the  Toxic 
Hazards  Research  Unit  (THRU).  The  THRU  is  the  contractor-operated  effort  of  the  Toxicology  Division, 
Occupational  and  Environmental  Health  Directorate,  Armstrong  Laboratory,  located  at  Wright-Patterson 
Air  Force  Base,  Ohio,  and  is  operated  under  Department  of  the  Air  Force  Contract  No.  F33615-90-0532. 
Allen  Vinegar,  Ph.D.,  Manager  of  the  THRU’s  Biological  Simulation  Section,  coordinated  this  effort. 
Lt  Col  Terry  A.  Childress  served  as  the  Contract  Technical  Monitor. 

The  research  described  herein  was  conducted  in  August  1994,  and  describes  factors  relevant  to  modeling 
the  uptake  of  chemical  toxicants  by  the  lung  during  brief  exposures  where  unsteady  state  events  prevail. 
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IMPROVEMENTS  IN  MODELING  OF  PULMONARY  UPTAKE  OF  TOXICANTS 


Richard  Collins 

ABSTRACT 


The  primary  objective  of  this  technical  report  is  to  provide  a  rational  foimdation  for  quantitatively 
evaluating  die  importance  of  unsteady  versus  steady  or  quasi-steady  events  during  uptake  of  chemical  toxicants 
via  the  human  respiratory  tract. 

A  secondary  objective  is  to  review  critically  the  current  approaches  in  physiologically  based 
pharmacokinetic  (PBPIQ  modeling  of  die  lung  widi  a  view  to  generalizing  the  model  descriptions  to  encompass 
a  broader  range  of  exposure  and  physiological  conditions. 

In  the  following  sections,  die  secondary  objective  is  addressed  first,  with  a  succinct  review  of  the 
transport  mechanisms  diat  can  operate  widiin  the  respiratory  tract  when  exposed  to  a  concentration  of  chemical 
toxicants  in  die  form  of  vapors,  aerosols,  and  particulate  matter.  Functional  relationships  are  reviewed  within 
die  context  of  a  voluminous  body  of  international  literature  bearing  on  this  and  related  subjects  of  importance. 

Computational  models  are  proposed  to  provide  an  essential  complementary  research  tool  with  not  only 
simulative  but  also  predictive  capabilities. 

Drawing  iqion  diis  considerable  background  of  established  knowledge,  one  can  proceed  more  readily 
widi  a  detailed  study  of  die  very  important  question  posed  as  the  primary  objective;  namely,  die  establishment 
of  an  order-of-magnitude  analysis  resulting  in  a  general  classification  scheme  that  identifies  three  major  flow 
regimes,  distinguished  on  the  basis  of  whether  the  flow  is  dominated  by  unsteadiness,  viscous  effects,  or  the 
effects  of  convective  acceleration. 

The  report  concludes  widi  some  important  recommendations  regarding  the  extrapolation  of  test  results 
performed  on  small  laboratory  animals  to  the  human  context. 

1.  INTRODUCTION 

Physiologically-based  pharmacokinetic  (PBPK)  modeling  of  die  pulmonary  uptake  of  vapors,  aerosols, 
and  particulates  draws  on  several  interrelated  scientific  disciplines.  These  include  biofluid  dynamics,  tissue 
mechanics,  biochemistry,  biophysics,  physiology,  and  mathematics,  culminating  in  the  mathematical 
formulation  and  computational  solutions  of  the  equations  of  mass  continuity,  momentum,  and  energy  as 
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embodied  in  nonlinear  simultaneous  systems  of  partial  differential  equations  with  associated  initial  and 
boundary  conditions. 

Multidisciplinary  teamwork  is  essential  for  such  detailed  investigations,  both  experimental  and 
tiieoretical,  into  the  underlying  mechanisms  governing  such  complex  human  biological  systems.  In  practice, 
as  much  e?q)erimentation  as  possible  is  performed  upon  animals  as  a  surrogate  for  humans.  Prior  to  initiating 
each  new  experimental  program,  it  is  useful  and  highly  advisable  to  formulate  a  series  of  models  of 
progressively  increasing  mathematical  complexity.  This  serves  a  dual  function  by  providing:  (a)  a  basis  for 
the  rational  design  of  an  experimental  progr  in  which  all  the  governing  biophysical  and  physiological 
parameters  and  variables  may  be  evaluated  to  prescribed  levels  of  accuracy,  and  (b)  a  framework  into  which 
tiiese  data  may  be  incorporated  and  interpreted.  Quantitative  comparisons  of  model  predictions  with  observed 
laboratory  measurements  serve  to  corroborate  the  realistic  fouiulation  of  successful  models,  permitting  the 
latter  to  be  usefully  employed  to  predict  responses  to  new  situations  in  lieu  of  conducting  further  costly 
laboratory  testing. 

In  die  event  tiiat  the  quantitative  comparisons  between  model  and  experiment  are  not  satisfactory,  the 
two  may  require  revision  through  an  interactive  process  until  the  predictions  of  the  model  reflect  the  reality 
of  the  (Aserved  phenomena  under  investigation.  This  process  is  universal  to  all  branches  of  scientific  research 
and  has  recently  been  summarized  in  the  context  of  PBPK  modeling  in  a  chapter  by  Leung  (1993). 

The  first  step  in  the  formulation  of  a  realistic  model  is  the  development  of  a  working  hypothesis  for 
die  presumed  mechanisms  of,  in  the  iM*esent  case,  uptake  of  toxicants  by  the  lungs  in  the  forms  of:  (a)  vapors, 
(b)  aerosols,  and  (c)  particulate  matter.  At  the  same  time,  based  on  logical  principles  of  biophysics  and 
physiology,  one  may  identify  the  primary  and  eventually  the  secondary  and  tertiary  parameters  and  variables 
diat  need  to  be  taken  into  account  in  such  a  model.  The  process  may  be  accomplished  in  successive  phases  of 
increasing  con^lexity,  gaining  valuable  experience  upon  the  completion  of  each  phase  of  the  overall  program. 
Alternatively,  one  may  enter  the  project  at  a  more  advanced  phase  at  the  outset  by  building  upon  and  by 
improving  the  elements  of  known  successful  models  by  either  refining  the  assumptions  made  in  those  earlier 
models  and/or  by  generalizing  or  adapting  them  to  the  new  situations  of  interest. 

With  this  general  framework  in  mind,  we  will  summarize  in  the  following  sections  the  general  nature 
of  PBPK  models  of  the  lung  previously  developed  for  purposes  of  predicting  pulmonary  uptake  imder  (quasi-) 
steady  state  conditions.  We  will  then  propose  refinements  of  flie  single  flow-limited  compartmental  model  to 
encompass  more  general  and  more  realistic  geometries,  possible  interactions  with  other  organ  systems,  and 
possible  ^neralizations  to  deal  with  acute,  short-term  exposures  to  toxicants  of  the  three  forms  cited  above. 
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2  PRESENT  METHODS  OF  MODELING  THE  LUNG 

In  currently  used  PBPK  models,  the  uptake  function  of  the  lung  is  idealized  in  die  form  of  a  single 
flow-limited  compartment  widi  a  unidirectional  alveolar  ventilation.  Uptake  of  vapors  depends  on  the  alveolar 
ventilation,  cardiac  ou^ut,  and  partitioning  of  tiie  intake  chemical  between- blood  and  air.  The  model  currently 
utilized  was  developed  using  the  general  structure  of  Ramsey  and  Anderson  (1984)  in  order  to  provide 
quantitative  descriptions  of  the  uptake,  distribution,  metabolism  and  elimination  of  xenobiotic  substances. 

Typically,  the  mass  balance  of  chemical  Al  in  the  lung  at  any  instant  in  time  is  described  by  a 

differential  equation  (Vinegar  et  al.  1990  b)  of  the  form: 

d  A,  =  Qajv  (Cjfjfj- Cgxu)  QeL(^r~^A) 
dt  (air  via  lungs)  (blood  via  lungs) 

+  PA  /(K^  +  CJ 

(air  across  lung  tissue)  (loss  due  to  metabolism) 

In  addition.  Vinegar  et  al.  (1990  a)  report  measurements  of  minute  alveolar  ventilation  rates  QPT 
which  were  then  correlated  for  an  anatomical  dead  space  Vp  by  a  term 

V^(t)  =  0.3  QPT(t)^f„^^ 

to  yield  an  effective  value  QP  of  alveolar  ventilation  in  flie  form 

QP(t)  =  QPT(t)-V^ 

Cardiac  ou^ut  was  tiien  assumed  proportional  to  the  effective  alveolar  ventilation  QP  in  the  form 

QC(t)  =Kc  •  QP(t) 

It  can  be  readily  corroborated  from  the  above  that  die  Ixmg  appears  to  have  been  modeled  in  its 
simplest  form,  that  of  a  single  rigid-walled  conqjartment  with  a  dead  space  volume  proportional  to  30%  of  tiie 
minimum  alveolar  ventilation  rate.  It  has  been  aptly  pointed  out  by  Gearheart  et  al.  (1993)  that  indirect 
estimation  of  (metabolic)  parameters  from  "fits"  to  a  PBPK  gas  uptake  model  will  depend  for  flieir  success 
on  the  adequacy  of  tiie  model  structure.  A  true  test  of  the  structural  adequacy,  and  indeed  tiie  usefulness  of 
the  model  will  be  its  general  capability  of  predicting  responses  accurately  over  a  sufficiently  wide  range  of 
parameter  values,  surpassing  that  of  the  laboratory  experiments  themselves. 

The  question  which  is  now  posed  is  to  ascertain  the  applicability  of  this  single  compartmental  model 
to  a  class  of  more  generalized  pulmonary  flows  for  which  uptake  of  toxicants  may  not  yet  have  reached  a 
"steady  state." 

Such  cases  would  relate  to  sudden  and  acute  pulmonary  exposure  to  a  toxicant  over  a  short  duration, 
or,  alternatively,  to  abrupt  cessation  of  and  subsequent  recovery  from  such  exposure.  In  either  case,  the 
pulmonary  respiration  and  circulation  undergo  a  substantial  transition  phase,  the  details  of  which  are  needed 
to  assess  quantitatively  tiie  associated  risks. 
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In  the  following  section,  we  will  return  to  first  principles,  examining  in  some  detail  flie  expected 
response  of  the  lung  to  highly  dynamic  changes  in  a  number  of  key  variables,  as  a  means  of  suggesting  that 
it  may  be  of  use  in  such  situations  to  envisage  a  modified  model  structure  which  more  fiilly  accounts  for  the 
physiological  interactions  which  will  come  into  play  in  this  more  complex  transitory  phase. 

Initially^  fte  prqwsed  improvements  in  the  PBPK  model  will  be  related  to  the  uptake  of  vapors,  and 
extended  only  later  to  tlte  more  complicated  uptake  of  liquid  aerosols  and  solid  particulate  matter. 

3.  GENERAL  OVERVIEW  OF  DYNAMIC  PULMONARY  VENTILATION  AND  VAPOR 

UPTAKE 

Ccmventionally,  die  lung  is  divided  into  an  alveolar  conqiartment  of  uniform  gas  concentration  wherein 
gas  exdiange  widi  die  blood  takes  place  by  diffusion,  and  a  dead  space  Vjj  in  which  no  gas  exchange  occurs. 

However,  ventilation  in  the  physiological  range  involves  mass  tranqiort  by  both  convection  and 
diffusion.  It  is  widely  accepted  tiiat  (bulk  flow)  convection  dominates  in  the  larger  diameter  airways  where 
flow  velocities  are  larger;  while  diffusion  is  more  significant  at  die  level  of  die  terminal  bronchioles  and 
alveoli,  where  vapor  transit  times  are  longer.  There,  mass  transport  develops  primarily  due  to  concentration 
gradients. 

PulmcHiary  airflow  profiles  in  the  largest  airways  are  flatter  dian  die  parabolic  profiles  characteristic 
of  fully  develqped  laminar  Poiseuille  flow.  Such  flows  are  likely  to  become  turbulent  as  they  pass  through 
successive  bifurcations  along  dto  tradieobrondiial  tree,  engendering  diaotic  eddy  currents  which  enhance 
mixing  befOTe  developing  a  more  laminar  character  as  the  flow  progresses  further  downstream.  This  process, 
which  is  well  known  in  classical  fluid  mechanics,  has  been  described  in  its  pulmonary  context  by  Hamilton 
(1986). 

The  beating  of  the  heart  further  contributes  to  airway  (cardiogenic)  mixing  as  the  medianical 
oscillations  of  die  heart  produce  vibration  of  the  walls  of  the  airways  (West  et  al.  1961).  Ihe  degree  of  mixing 
which  results  may  be  taken  as  largely  independent  of  heart  rate  (Dreschler  et  al.  1984).  Indeed,  such  cardiac 
oscillations  can  theoretically  enhance  gas  mixing  in  the  airway  of  the  dog  by  1000%  (Slutsky  1981)  dirough 
the  coupled  effects  of  molecular  diffusion  and  augmented  convective  transport. 

In  the  smaller  airways,  where  the  flow  is  laminar,  Taylor  di^rsion  promotes  enhanced  radial 
diffusive  mWing  in  adjacent  annular  rings  as  a  result  of  local  concentration  gradients,  in  addition  to  axial 
mixing  along  the  length  of  the  airway.  This  tends  to  produce  a  bi-directional  flow:  the  components  witii  the 
hi^r  alveolar  concentration  are  tranqiorted  outward  along  die  walls  of  the  airway,  while  die  components  widi 
the  higher  mouth  concentration  appear  to  move  toward  the  alveoli  along  the  center  core  of  die  airway. 
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The  process  is  qrclic  and  hi^y  (fynamic.  Following  an  inspiration,  the  “tongue”  of  the  Oj  pulse  from 
die  moufli  extends  down  die  airway,  diffusing  out  to  the  wall,  while  at  the  same  time,  CO2  diffuses  in  toward 
die  center  of  die  airway.  The  latter  C02-rich  air  is  then  carried  out  toward  the  mouth  during  expiration,  while 
O2  diffuses  into  the  tongue  of  die  pulse  as  COj  diffuses  outward.  The  process  is  repeated  at  each  new  cycle, 
and  is  cumulative  so  that  totol  mixing  is  augmented,  particularly  at  higgler  frequencies,  even  for  small  tidal 
volumes  less  than  the  dead  space. 

It  has  been  proposed  by  Lehr  (1980)  tiiat  die  different  regions  of  the  lung  do  not  inflate  uniformly,  but 
may  at  any  time  be  up  to  180°  out  of  phase  with  one  anodier.  Indeed,  stroboscopic  photography  of  excised  dog 
lungs  indicates  a  “Pendelluft”  movement  with  significantly  different  time  constants  (resistance  x  compliance) 
for  different  individual  airway  units,  engendering  retrograde  flows.  Very  similar  predictions  came  out  of  die 
computations  of  Collins  et  al.  (1979)  for  blood  flow  in  die  lung  as  a  direct  consequence  of  significant  measured 
differences  in  die  anatomical  structures  and  wall  compliances  of  die  pulmonary  vasculature.  They  successfully 
mathematically  modeled  the  flow  profiles  within  die  five  pulmonary  interconnected  lobes  of  the  pulmonary 
circulation. 

The  asymmetry  of  the  airway  architecture  and  material  properties  can  Induce  gas  exchange  in  regions 
whose  anatomical  dead  space  Vd  is  less  dian  the  global  dead  space.  In  such  regions  with  tidal  volumes  Vt  < 
Vd,  adequate  alveolar  ventilation  could  obtain  widiout  resorting  to  diffusive  transport. 

Steady-state  exdiange  is  proportional  to  tidal  volume  Vt  minus  dead  space  volume  Vq  and  no  steady 
state  gas  exchange  occurs  for  Vt  ^  V^.  According  to  Mitzner  (1986),  diis  conclusion  is  based  on  the 
assumption  diat  gases  traverse  die  dead  space  widi  a  blunt  flow  profile  (velocity  constant  over  the  airway  cross- 
section)  during  die  breathing  (tycle.  Under  such  conditions,  a  small  tidal  volume  Vt  could  not  traverse  a  large 
dead  space  volume  Vq.  However,  die  velocity  profile  is  not  blunt  or  flat,  but  rather  tongue-shaped,  as 
discussed  earlier.  Thus,  in  diis  context,  some  a^reciable  alveolar  ventilations  and  gas  exchange  can  occur  for 
Vt  <  Vd,  a  result  which  could  not  be  modeled  readily  using  a  simple  one-  or  two-compartment  idealization 
of  a  steady-state  pulmonary  airflow. 

In  fact,  it  would  appear  diat  this  bi-directional  axial  dispersion,  or  spreading  of  die  gas  front  into  a 
tongue  along  die  tube  (airway)  axis  during  both  inspiration  and  expiration,  contributes  significandy  to  die  gas 
mixing  process.  A  corresponding  model  then  should  include  branching  networks  of  idealized  airway  tube 
segments,  with  varying  airway  padi  lengdis  and  corresponding  transit  times.  Scbroter  and  Sudlow  (1969)  have 
measured  significandy  different  flow  velocity  profiles  during  inspiration  and  expiration  in  a  physical  branched 
airway  model. 
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Accordingly,  an  improved  global  model  is  proposed  for  vapor  uptake  reflecting  this  inherent 
asymmetry  over  the  the  two  halves  of  the  breathing  cycle,  as  flow  streamlines  cannot  be  precisely  retraced 
between  inspiration  and  expiration.  TTie  true  anatomical  dimensions,  branching  characteristics  and  physical 
properties  of  the  distensible  airways  may  be  used  to  characterize  the  xmderlying  mechanisms  responsible  for 
the  overall  axial  dispersion  of  the  flow  profiles.  From  these,  realistic  gas  mixing  rates  can  be  computed  for 
a  variety  of  unsteady  pulmonary  flow  regimes. 

4.  DIFFUSION  AND  SETTLING  OF  AEKOSOL  PARTICLES  IN  THE  PULMONARY 
AIRWAYS 

According  to  Hext  et  al.  (1993),  inhalation  represents  one  of  the  major  routes  by  which  humans  can 
be  exposed  to  foreign  materials.  These  may  then  be  readUy  absorbed  or  may  react  directly  with  adjacent  tissues 
within  the  body.  The  successful  experimentalist  requires  a  quantitative  understanding  of  the  anatomy  of  die 
respiratory  tract,  die  mediainans  of  deposition  of  liqitid  (aerosol)  or  solid  (particulate)  matter  and  the  species 

differences  between  die  experimental  animal  and  humans. 

It  is  quite  common  for  irritant  toxic  aerosols  and  particulate  matter  to  influence  considerably  the 
breathing  patterns  of  test  animals,  possibly  reducing  the  re^iratory  rate  1^  80%  This  reduction  can  be 
concentratirai-related  and  response  relationships  can  be  plotted,  for  example,  as  a  function  of  reflex  response 
and/or  bronchoconstriction.  For  liquid  aerosols,  solubility  in  humidified  saturated  atmo^ieres  must  be 
accounted  for. 

4.1  Transport  Medianisins 

Dqx^ticm  of  particles  widiin  tte  respiratory  tract  is  governed  Ity  five  principal  mechamsms:  inertial 
in^iaction,  gravitatiraial  sedimentatk®.  Brownian  diffurirai  and  interception  and  electrostatic  precipitation.  Hie 
particle  size  and  density  relative  to  air  are  important  in  assessing  die  relative  roles  of  diese  mechamsms  in 
various  regions  of  the  respiratory  tract. 

Experimental  results  indicate  diat,  for  a  single-breadi  delivery,  the  fraction  of  inhaled  aerosol 
recovered  closely  follows  an  exponential  function  of  breath-holding  time  for  particle  diameters  ranging  from 
0.1/1  to  0.8/1.  Multi-expraiential  recovery  curves  are  required  for  larger  particles  (Goldberg  et  al.  1981).  Axial 
di^rsion  effects  were  incorporated  by  Yu  et  al.  (1979)  into  their  model  for  the  analysis  of  loss  of  aerosol  at 
different  depths  widiin  die  respiratory  tract  during  inhalation  and  exhalation. 


6 


Larger  particles  of  density  different  than  air  will  not  follow  the  flow  streamlines  of  the  inhaled 
pulmonary  airflow.  They  will  tend  to  deposit  in  the  upper  respiratory  bronchi,  where  they  can  be  moved 
mouthward  into  the  esophagus  by  mucociliary  transport. 

4.2  Mucociliary  Clearance 

It  has  been  observed  fliat  peak  e3q)iratory  flow  rates  are  systematically  greater  than  corresponding  peak 
inspiratory  rates  (Chang  1986).  These  different  flow  rates  produce  biased  air-mucus  shear  forces  favoring 
moutiiward  movement  of  mucus  or  sputum.  But  an  important  role  may  also  be  played  by  vagal  stimulation  of 
ciliary  beating  via  efferent  receptors.  Such  clearance  mechanisms  appear  to  operate  during  coughing. 
Enhancement  of  tracheal  mucus  clearance  rate  (TMCR)  is  most  pronounced  for  breathing  frequencies  in  the 
range  of  1 1  -  15  Hz  in  anestiietized  dogs,  reaching  peak  values  of  340%  of  control  at  13  Hz  (King  et  al.  1984) . 
The  velocity  of  mucus  movement  is  proportional  to 

dVrfUR-1) 

where  6  is  flie  titickness  of  the  mucus  layer,  Vj  is  the  tidal  volume,  f  is  frequency,  and  R  is  die  ratio  of  peak 
expiratory  to  peak  inspiratory  flow  rate  (Chang  et  al.  1987). 

5.  FORMULATION  OF  A  PROPOSED  COMPUTATIONAL  MODEL  FOR  HUMAN 
PULMONARY  DYNAMICS 

A  con^lete  branching  model  of  tiie  respiratory  system  is  proposed  on  the  basis  of  detailed  anatomical 
measurements  (Fredberg  et  al.  1978).  A  quasi  one-dimensional  formulation  for  the  coupled  fluid/wall 
interaction  can  be  resolved  numerically  using  a  modified  two-step  Lax-Wendroff  finite-difference  technique. 
This  model  can  subsequently  be  used  to  estimate  the  regional  deposition  of  fine  aerosol  particles  in  the  human 
lung. 

Preliminary  theoretical  models  of  particle  deposition  coupled  with  experimental  studies  appear  to 
corroborate  the  overall  importance  of  particle  size  (Hext  et  al.  1993).  Particle-laden  air,  inhaled  through  the 
nose,  finds  its  way  through  to  the  ends  of  the  trachea,  where  the  branching  of  the  airways  begins;  extending 
eventually  down  through  the  bronchi  and  bronchioles  to  the  alveolar  ducts  and  sacs.  In  the  Landahl  (1950) 
model  of  the  respiratory  tract,  eleven  primary  regions  have  been  identified,  along  with  their  average 
anatomical  dimensions  and  corresponding  air  velocities.  It  will  remain  to  associate  with  these  tiie 
corresponding  mg^hanical  properties  of  wall  elasticity  or  compliance  for  use  in  tiie  mathematical  model  to  be 
formulated. 
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As  the  particles  entrained  in  the  airflow  do  not  have  die  same  density  as  the  surrounding  air,  their 
trajectories  will  not  follow  die  aerodynamic  streamlines.  Particularly  at  bifurcations  and  flow  dividers,  where 
the  curved  streamlines  predispose  to  the  development  of  secondary  flows,  the  centrifugal  forces  on  such 
particles  will  push  them  laterally  out  toward  the  walls  where  they  may  deposit  by  impact  and  be  captured. 

As  the  airways  continue  their  branching  with  successive  generations,  the  cross-sectional  area  of  the 
respiratory  network  increases  dramatically.  Air  velocities  decrease  from  150  cm/sec  at  die  trachea  to  about 
2  cm/sec  at  the  terminal  brondiioles.  It  is  e^cially  in  this  region  of  slowly  moving  flow  that  one  would  expect 
an  increased  influence  on  particle  deposition  due  to  the  action  of  sedimentation  and  Brownian  diffusion. 

To  a  lesser  extent,  electrical  charges  associated  with  particles  dispersed  into  the  air  by  grinding  or 
combustion  may  cause  particle  aggregation  and  a  concomitant  alteration  in  the  effective  particle  size 
distribution,  with  their  consequences  on  particle  deposition  patterns. 

Thermal  forces  are  likely  even  less  important.  Althougjh  the  temperature  in  die  respiratory  tract  may 
rise  sli^dy  dining  inhalation,  these  minimal  temperature  differences  are  thought  to  have  only  a  minor  effect 
upon  particle  deposition  due  to  Brownian  motion  (Morrow  1960). 

5.1  Mathanatical  formulatuni 

The  first  task  is  to  determine  die  flow  patterns  of  air  in  the  re^iratory  network.  The  fluid  equations 
of  mass  continuity  and  momentum  must  be  solved  in  conjunction  with  a  realistic  “law  of  the  wall”  which 
relates  local  cross-section  of  the  airways  to  instantaneous  transmural  pressures  as  a  function  of  local  wall 
properties  (modulus  of  elasticity).  Special  provisions  can  be  made  for  airway  collapse  under  qioradic 
conditions  of  negative  transmural  pressure. 

The  governing  equations  are  solved  numerically  in  a  manner  similar  to  that  used  by  Collins  et  al. 
(1979),  for  die  pulmonary  circulation.  The  true  physiological  variations  in  pressure  and  flow  rate  as  recorded 
during  respiration  will  serve  as  boundary  conditions  to  the  model.  Once  die  flow  field  has  been  computed  in 
this  manner,  one  may  calculate  the  transport  of  particles  of  different  sizes  and  densities  which  are  entrained 
in  this  known  flow.  This  proposed  decomposition  of  the  problem  tacitiy  implies: 

a)  die  basic  “clean  air”  flow  will  not  be  substantially  altered  by  die  intake  of  dispersed  particles, 
and 

b)  a  lack  of  interaction  between  particles  of  differing  sizes  as  diey  deposit  onto  the  walls  of  the 
respiratory  tract.  TWs  allows  die  deposition  distributions  to  be  summed  for  given  particle  size 
distributions. 
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The  alternative  to  fliis  more  simplified  approach  would  be  a  fully  coupled  two-phase  treatment  of  the 
particle-air  equations,  including  particle-particle  interactions  in  a  non  stationary  two-dimensional  flow  with 
rotational  symmetry.  The  zones  of  possible  momentary  turbulence  and  flow  re-circulation  can  be  evaluated 
in  terms  of  corresponding  losses  in  stagnation  pressure  due  to  viscous  dissipation  for  given  bifurcation  angles 
and  parent-daughter  area  ratios.  Nondichotomous  branching  may  impose  additional  difficulties,  leading  to  the 
necessity  of  discretizing  groups  of  similar  airway  generations  into  equivalent  simplified  networks. 

Such  an  approach  should  constitute  a  substantial  improvement  over  the  analysis  of  Taulbee  and  Yu 
(1975),  wherein  the  airway  geometry  was  presumed  to  be  boimded  by  rigid  walls.  Lanunar  flow  was 
considered  exclusively,  and  axial  diffusion  and  particle  deposition  in  the  mouth  and  trachea  were  neglected. 
The  effects  of  airway  bifurcations  on  viscous  dissipation  and  total  pressure  losses  were  not  dealt  with  in  that 
formulation. 

6.  ORDER-OF-MAGNTTUDE  ANALYSIS  OF  PULMONARY  FLOW  PARAMETERS: 

RELATIVE  IMPORTANCE  OF  UNSTEADINESS 

Pulmonary  airflow  and  the  response  of  the  respiratory  system  to  die  uptake  of  toxicants  constitute  a 
hi^y  complex  biological  system  which  is  extremely  difficult  to  analyze.  Within  an  intricate  and  truly  three- 
dimensional  network  of  branching  distensible  conduits,  oscillatory  flow  phenomena  occur  over  a  wide 
spectrum  of  interacting  frequencies.  It  is  dierefore  not  surprising  that  applied  bio-madiematicians  seek  to 
capture  the  primary  flow  mechanisms  within  a  simpler  analogous  system  which  is  more  tractable  to 
computational  solution  and  which,  nonetheless,  retains  the  salient  features  of  the  observed  physiological 
response  (Jan  et  al.  1989). 

In  particular,  it  is  die  presence  of  time  derivatives  (unsteadiness)  in  the  governing  mathematical 
equations  diat  gready  complicates  die  solution.  If  one  could,  for  example,  determine  flow  regions  over  which 
die  unsteacfy  terms  played  a  lesser  role  than  do  odier  terms  in  the  same  transport  equations,  then  the  flow  could 
be  considered  as  quasi-steady,  permitting  considerable  simplification  of  the  mathematical  solution.  By  quasi- 
steady,  one  means  that  notwithstanding  the  time-varying  flow  rates,  the  detailed  characteristics  of  the  flow  at 
any  instant  (snapshot)  in  time  “adjust”  themselves  immediately  to  temporal  changes  in  the  external  flow 
conditions.  Flow  fields  dien  would  be  identical  to  those  that  would  obtain  in  a  truly  steady  flow  at  the  same 
flow  rate  at  that  instant  of  time.  In  odier  words,  no  “inertial”  influence  or  time  lag  due  to  the  unsteadiness  of 
the  flow  would  be  discernible. 

Pedley  et  al.  (1986)  have  shown  diat  for  a  lung  idealized  as  a  branching  network  of  rigid,  self-similar 
bifurcations,  the  purely  oscillatory  flow  of  an  incompressible  fluid  may  be  characterized  by  just  two 
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independent  dimensionless  groups;  namely,  a  dimensionless  frequency  and  a  stroke  length  parameter  L/a, 
where  a  is  die  airway  radius,  v  die  kinematic  viscosity  and  w  the  oscillatory  frequency.  L  is  die  average  axial 
Higpiarfttnftnt  of  a  fluid  particle.  Almost  all  other  dimensionless  groups  found  in  the  literature  can  be  formed 
from  combinations  of  these  two.  They  have  been  summarized  in  Table  1  appearing  in  the  Appendix  (Jan  et 
al.  1989).  For  the  human  airways,  with  normal  breathing  at  relatively  low  frequencies,  the  dimensionless 
frequency  ranges  between  0.01  <  a  <  3,  with  amplitudes  between  100  <  L/a  <  2000. 

Over  recent  years,  a  parallel  bo<fy  of  relevant  literature  has  developed  in  response  to  interest  in  high- 
frequency  ventilation,  characterized  by  0.2  <  a  <  20  and  20  <  L/a  <  500.  Indeed,  almost  80  years  ago, 
Hendersen  et  al.  (1915)  demonstrated  that  fresh  air  could  reach  the  other  end  of  the  respiratory  tract  with  a 
smaller  flushing  volume  than  the  anatomical  dead  space.  This  high-frequency  regime  is  of  interest  because 
stroke  volumes  can  be  reduced  by  a  factor  of  4  or  5  whUe  die  breathing  frequency  rises  by  factors  of  7  to  20. 
Study  of  this  regime  was  motivated  by  the  need  to  provide  assisted  ventilation  in  die  treatment  of  various 
respiratory  pathologies  characterized  by  reduced  tidal  volumes  and  low  peak  inspiratory  airway  pressures.  It 
is  not  unexpected  that  investigations  of  high-frequency  ventilation  have  spurred  on  developments  leading  to 
a  better  understanding  of  die  physiological  range  of  respiratory  parameters  as  well. 

Pedley  et  al.  (1986)  have  reviewed  die  earlier  literature  on  pulmonary  fluid  dynamics  in  which  steady 
and  quaa-steady  in^iratctfy  flows  were  studied.  Without  entering  into  the  full  technical  detaUs  in  this  succinct 
report  concerning  the  classical  theory  of  unsteady  flows  of  a  viscous  compressible  fluid,  the  main  physical 
ideas  will  be  presented  for  readers  not  having  a  specialized  background  in  die  field. 

As  flow  develof®  widiin  a  cylindrical  tube  sudi  as  a  pulmonary  bronchus,  the  air  velocity  varies  across 
die  tube,  fr<Mn  zero  at  die  walls  (to  which  the  air  particles  adhere  due  to  viscosity)  to  a  maximum  value  at  die 
tube  axis.  The  more  slowly  moving  fluid  layer  adjacent  to  die  tube  walls  forms  a  boundary  layer,  which  may 
be  diin  relative  to  the  tube  diameter,  if  the  Reynolds  number  Re  (ratio  of  inertial  to  viscous  forces)  is  high,  or 
may  fill  the  tube  progressively  as  in  fully-developed  Poiseuille  flow  characterized  by  smaller  Reynolds 
numbers. 

Beneath  diis  steady  (Blaaus)  viscous  boundary  layer  may  lie  a  thinner  unsteady  (Stokes)  viscous  layer. 
These  boundary  layers  are  of  qiecial  interest  at  bifurcations  where  the  flow  streamlines  curve  in  passing  from 
a  parent  to  a  daughter  conduit  during  inspiration  as  well  as  the  reverse  during  expiration.  At  such  bifurcations, 
most  of  the  pressure  losses  and  gas  mixing  occur. 
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It  is  suggested  that  a  parameter  characterizing  the  intrinsic  degree  of  unsteadiness  of  a  flow  may  be 
based  on  the  ratio 


e  ~L^/L 

where  is  die  lengdi  of  a  generation  of  conduits  and  L  is  the  average  axial  displacement  of  a  fluid  particle, 
and  is  equal  to  the  stroke  volume  divided  by  the  local  cross-sectional  area.  Such  a  parameter  €  or  its  variant 
e*  may  be  interpreted  physically  as  the  ratio  of  either: 

[€*]  (i)  unsteady  (local)  acceleration  to  the  convective  acceleration,  or 

[e/4]  (ii)  patii  length  of  a  generation  of  conduits  to  tiie  fluid  displacement  over  an  oscillatory 

half-cycle  (inspiration/expiration),  or 

[eVi]  (iii)  time  required  for  tiie  thickness  of  the  unsteady  Stokes  viscous  boundary  layer  to 
reach  the  thickness  of  the  steady  Blasius  viscous  layer  associated  commonly  witii  Poiseuille 
flow. 

The  nonlinear  differential  equations  governing  fluid  motions  are  given  in  vector  notation  in  the 
Appendix  and  will  not  be  repeated  here.  They  represent;  (a)  conservation  of  mass  of  an  incompressible  fluid 
and  (b)  conservation  of  momentum.  The  latter  is  simply  a  statement  of  Newton's  law;  force  =  mass  x 
acceleration.  The  forces  are  viscous  friction  (between  adjacent  particles  and  between  air  particles  and  the 
conduit  walls)  and  centrifugal  (forces  transverse  to  curved  streamlines)  in  addition  to  tiie  pressure  gradients 
present  in  the  flow.  The  acceleration  has  both  tenporal  (measured  at  a  fixed  position)  and  convective 
(measured  at  a  fixed  time)  components,  as  indicated  in  Table  1  of  tiie  Appendix. 

By  comparing  like  magnitudes  of  these  terms  in  a  classical  manner,  one  is  able  to  establish  the 
existence  of  at  least  three  flow  regimes  as  depicted  in  Fig.  1  (Appendix).  The  details  given  in  tiie  Appendix 
will  be  self-explanatory  for  a  reader  having  an  introductory  backgroimd  in  fluid  mechanics. 

Jan  et  al.  (1989)  have  carried  out  detailed  flow  experiments  on  a  plexiglass  model  of  a  pulmonary 
bifurcation,  measuring  the  full  3-dimensional  unsteady  flow  field  for  a  wide  variety  of  flow  conditions.  The 
results  will  be  (simplistically)  summarized  below  for  tiie  non-specialized  reader. 

To  ascertain  wheflier  a  flow  is  quaa-steady,  one  must  consider  how  quickly  that  flow  will  “re-adjust” 
or  become  compatible  witii  a  sudden  time-diange  in  the  pressure  gradient.  In  tiie  absence  of  swirling  secondary 
flows  witii  curved  tubes  (centrifugal  forces),  the  flow  can  only  readjust  to  a  change  in  the  pressure 

gradient  by  diffusing  vorticity  (or  equivalently,  momentum)  through  viscous  forces  over  the  cross-section  of 
tiie  airway.  The  time  required  is  appronmately  a% .  Witii  tiie  addition  of  the  secondary  flows  described  above 
for  curved  streamlines,  tiie  flow  can  adjust  more  quickly  because  of  the  associated  convective  mixing.  This 
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speeds  xq)  ibe  iaro<»ss  to  a  time  of  the  order  of  aA^s-  Here,  a  is  the  airway  radius,  v  is  the  kiiiematic  viscosity 
and  Vs  is  flie  cross-flow  velocity.  Since  a^/v  >  aA^s»  the  flow  is  certainly  quasi-steady  whenever  the  larger 
possible  value  of  die  readjustn^nt  tune 

«Vv  <  1/m 


or  alternatively,  when 

ec  • - <  1 

V 

This  implies  diat  flow  Can  be  comidered  to  be  quasi-steady  in  parts  of  convective  flow  regime  m 
jttctured  in  Fig.  1  (y^jpendix).  Implications  for  die  quasi-steady  uptake  of  gases  in  this  region  HI  will  now  be 
discussed  very  briefly. 

As  flow  reverses  during  the  tranaticm  fircmi  inspiration  to  expiration  (and  vice  versa),  significant  fluid 
miTirtg  occurs  within  the  redrculating  patterns  of  vorticity  which  form  momentarily  near  the  walls  of  die  flow 
divider.  Ideally,  maTimnm  lateral  mixing  over  the  complete  tube  cross-section  will  follow  immediately  upon 
die  reversal  in  direction  of  an  axial  peaked  "tongue"  of  convected  fluid.  An  effective  axial  diffusivity  can  be 
formulated  for  such  a  tranqiort  process  that  varies  as  the  square  of  the  tidal  volume  times  the  oscillation 
frequency  (Permutt  et  al.  1985). 

6.1  PULMONARY  MODELING  IN  REGIONS  I,  n  and  m 

It  is  suggested  that  flows  in  die  unsteady  zoi^  I  (Fig.  1)  for  which  stroke  lengdis  are  short  and 
secondary  motitms  small,  be  modeled  according  to  the  now  classical  linearized  theory  of  Womersley  (1955). 
One  can  assume  PoiseuUle-type  flow  everywhere  in  die  viscous  zone  n  where  secondary  flows  are  negligible 
and  die  Reynolds  number  Re  <30  (indicating  very  litde  inertia). 

The  ctMivective  zone  in  may  be  fiirdier  subdivi^  rcaighly  into  a  very  complex  turbulent  region  whidi 
has  not  yet  been  nxxleled.  In  this  ^bregion,  unsteady  eflects,  though  mit  dominant,  may  nonetheless  flare  up 
during  brief  periods  of  flow  acceleration.  They  are,  however,  much  weaker  at  peak  flow  or  during 
(teceleiative  (diases.  Outade  diat  aibre^on,  inqmtoiy  flow  could  be  modeled  as  an  entrance-length  transition 
ov^  ^ihidi  boundary  layers  gradually  build  up.  Expiratory  flow  (daughter-to-parent  conduits)  remains  to  be 
modeled. 
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6.2  UNSTEADY  FLOW  REGIME 

Unsteady  transverse  secondary  flows  contribute  significantly  to  the  mixing  of  chemical  toxicants  with 
die  respiratory  airflow,  particularly  at  the  sites  of  tracheobronchial  bifurcations  at  the  end  and  onset  of  each 
inhalation  and  exhalation.  This  mixing  process  is  itself  critically  determined  by  a  delicate  interplay  of  tiie 
centrifugal  forces  (U^/R)  on  the  fluid  or  solid  particles  and  the  transverse  components  of  the  temporal 
acceleration. 

For  unsteady  effects  to  dominate  tiiis  mixing  process,  one  requires  (see  Eq.  8,  9)  fliat  a^»  1.  This 
Womersley  number  will  now  be  evaluated  for  typical  airflows  in  flie  human  lung  at  normal  breathing 
frequencies  and  for  accelerated  (hi^  frequency)  ventilation. 

Using  an  average  value  of  the  kineniatic  viscosity  v  of  air  equal  to  1.5  x  10’*  m^/sec,  one  can  readily 
compute  the  Womersley  number 

0^  =  0^  (i)/v 

for  a  range  of  values  of  airway  radius  a.  Referring  to  flie  human  morphometric  model  (whole  lung)  of  Overton 
(1990)  (his  Appendix  B,  p.  317),  one  notes  a  wide  range  of  values  of  radius  a  from  1  cm  at  generation  1  to 
1/4  at  generation  7,  decreasing  down  to  0.015  at  tiie  terminal  bronchioles  of  generation  25.  Typical  values 
of  at  generation  g  are  presented  in  Table  2  below:  (normal  breathing  frequencies:  rat  =  75/niin;  human ' 
=  12/ntinute) 


TABLE  2. 

ESTIMATES  OF  THE  WOMERSLEY  NUMBER  AT  VARIOUS  GENERATIONS 
(BrtJMAN) 

Generation 

Airway  Radius 

Breathing  Frequencies 

Womersley  Number 

g 

a  (meters) 

(D  Hx 

Trachea 

1 

1.0  xlO’* 

10,  1,  0.2 

100,  10,  5,  2 

5 

3.3  xlO’* 

1,0.2 

0.67,  0.134 

6 

2.87  X  Iff* 

1,0.2 

0.55,  0.11 

7 

2.17x10’* 

1,0.2 

0.31,  0.062 

25 

3.0  xKr* 

1,  0.2 

6x10’*,  1.2xlff* 

Terminal  Bronchioles 

From  tiiese  preliminary  estimates,  it  is  clear  that  unsteadiness  (a*»  1)  in  human  respiration  at  normal 
breathing  frequencies  is  inherently  evident  in  the  first  few  generations  distal  to  the  trachea.  As  breathing 
frequency  increases  as  a  result  of  work,  exercise  or  neurogenic  factors,  etc.,  the  importance  of  unsteadiness 
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will  extend  down  to  the  hi^er  generation  airway  segments  as  well.  This  extension  is  readily  estimated  by 
noting  that  increases  linearly  with  breathing  frequency  u . 

In  contrast  to  die  human,  die  Womersley  number  for  the  rat  is  shown  in  Table  3  to  be  consistently 
much  less  than  unity  at  normal  breadiing  frequencies  for  all  generations  of  the  respiratory  tract,  including  the 
first  few  generations  distal  to  the  trachea. 


TABLE  3.  Estimates  of  the  Wommley  number  at  various  gMierations  (rat) 


Generation 

g 

Airway  Radius 
a  (met^) 

Breathing  Frequencies 
(i)  Hz 

Womersley  Number 

Of* 

**  g 

Trachea 

1 

1.7x10'’ 

1 

0.193 

5 

8.1X10-^ 

1 

4.39x10-’ 

6 

6.7x10"' 

1 

3.00x10-’ 

7 

ti.lxlO"* 

1 

2.49x10’ 

25 

Terminal  Bronchioles 

4.3x10-’ 

1 

1.23x10'’ 

This  difference  with  human  response  is  due  to  smaller  airway  diameters  for  the  rat.  These  diameters 
ai^iear  stpiared  in  die  exp'esson  for  die  Womersley  number  and  therefore  outweigh  the  approximately  5-fold 
increase  in  the  rat's  normal  breathing  frequency  (1  Hz)  relative  to  die  human  (0.2  Hz)  used  in  diese 
calculations. 

This  important  finding  would  indicate  that  respiratory  airflow  is  essentially  quasi-steady  throughout 
the  totality  of  the  respiratory  tract  of  the  rat,  whereas  this  is  not  so  for  the  human.  Accordmgly,  respiratory 
testing  of  toxicant  uptake  in  the  rat  may  not  always  display  tiie  same  qualitative  flow  regimes  observable  in 
die  human  upper  airways  (s^  Rg.  1).  Extrapolations  of  data  taken  on  die  rat  to  assess  human  response  should 
take  these  significant  differences  into  account  in  an  aj^ropriate  manitor. 

In  die  models  develc^ied  by  Ramsey  and  Andersen  (1984)  and  by  Overton  (1990),  the  anatomy  of  the 
pulmonary  airways  has  been  “lumped”  into  at  least  three  compartments  associated  with  die  upper  respiratory 
tract  (URT),  the  tracheobronchial  (TB)  region  and  the  pulmonary  region,  rei^ctively. 

Furdiermore,  good  agreement  widi  experimental  data  was  achieved  by  Ramsey  and  Andersen  (1984) 
only  after  fitting  die  model  simultaneously  to  data  from  four  experiments  in  order  to  determine  compatible  sets 
of  physiological  parameters.  Even  diou^  such  estimated  model  parameters  may  be  reasonable,  it  is  important 
to  verify  that  they  will  hold  also  for  widely  different  experimental  data  to  which  they  have  not  been  fitted.  The 
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true  value  of  a  model  lies  in  its  ability  to  predict  response  to  a  wide  range  of  operating  conditions  well  beyond 

those  envisaged  in  its  formulation  wifli  experiment. 

Generalized  models  can  certainly  be  improved  substantially  for  more  accurate  predictions  of  sudden 
exposure  and/or  cessation  of  exposure  to  chemical  toxicants  by  incorporating  a  realistic  unsteady  analysis  into 
fee  URT  and  possibly  the  TB  compartments  of  the  Overton  model.  In  addition,  realistic  geometry  reflecting 
tile  anatnmiral  details  of  the  respiratory  tract  should  be  incorporated  into  the  models  at  all  levels,  if  possible. 

Fmally,  basic  physiological  parameters  should  be  estimated  by  direct  measurement  wherever  possible* 
and  not  obtained  flirou^  fitting  of  model  results  to  experimental  data.  If  several  parameters  are  fitted 
simultaneously  in  such  a  process,  it  is  not  at  all  guaranteed  that  the  parameter  determination  will  be  unique! 

6.3  IMPLICATIONS  FOR  THE  EXTRAPOLATION  OF  TOXICOLOGICAL  ANIMAL  TEST 
RESULTS  TO  HUMANS 

From  tiiis  very  preliminary  survey  of  modeling  of  pulmonary  airflows,  it  becomes  abundantiy  clear 

that: 

1)  current  quasi-steady  PBPK  models  fi)r  estimation  of  the  uptake  of  chemical  toxicants  in  small 
animals,  typically  rats,  must  correspond  to  tiie  particular  flow  regime(s)  prevailing  in  flie 
animals, 

2)  tile  flow  regimes  applicable  to  human  pulmonary  iqitake  of  given  toxicants  must  exist  in  tiie 
animal  (ests  as  a  necessary,  if  not  sufficient,  condition  for  evaluating  quantitatively  the  same 
toxicant  transport  mechanisms  of  interest  to  human  uptake. 

rnmpiiance  with  these  essential  requirements  may,  in  some  cases,  require  careful  reconsideration  of 

the  basic  design  of  particular  animal  testing  programs. 

6.4  FUTURE  WORK 

It  would  be  very  wortiiwtule  to  formulate  a  detailed  respiratory  model  of  the  human  lung  based  upon 
a  branching  network  of  distenable  pulmonary  airways.  Species  equations  for  each  of  tiie  chemical  constituents 
can  then  be  incorporated  and  metabolic  rate  functions  assigned  locally.  Model  results  may  flien  be  compared 
quantitatively  with  those  of  tiie  previous  compartmental  models  as  a  standard  for  evaluation  of  accuracy  in 
ajcwssing  risk  due  to  uptake  of  chemical  toxicants  via  tiie  respiratory  tract. 
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APPENDIX 


NONSTATIONARY  FLOW  IN  THE  LUNG 


Order-of-Magnitude  Analysis  of  Dominant  Flow  Regimes 
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NONSTATIONARY  FLOW  IN  THE  LUNG 

In  order  to  interpret  the  qualitative  features  of  different  flow  regimes  which  may  develop  within  the 
respiratory  tract,  it  is  very  helpful  to  formulate  an  overall  classification  scheme  based  upon  an  order-of- 
magnitude  analysis  of  die  underlying  fluid  dynamic  equations.  These  equations,  representing  conservation  of 
mass  and  momentum,  are  constituted  by  a  system  of  nonlinear  partial  differential  equations  with  initial  and 
boundary  conditions.  For  tiie  purposes  of  this  preliminary  analysis,  these  will  be  expressed  below  for  axially 
symmetric  flow: 

<frv  F  =  0 

dV  “^1 

—  +  VgradV  «  —  grad p  +  v  divgradV 
dt  p 


where  V  is  the  velocity  vector  and  p  the  pressure  in  the  airflow.  ITie  second  vector  equation  represents  a 
simple  force  balance  between  temporal  acceleration,  convective  acceleration,  pressure  gradient  aiKi  viscous 
forces,  re^ctively.  For  axially  symmetric  flows,  they  reduce,  without  loss  of  order-of-magnitude  generality, 
to 

dx  ^ 


du  ^  ^  du  ^  du  1^ 

^  dx  dy  p  dx 


di>  du  I  dp  6^0 

dt  dx  ^  p  dy  Sy'^ 


where  V(/x,v)  denotes  the  axial  and  transverse  components  of  the  fluid  velocity  respectively,  in  the  x  and  y 
directions.  These  various  terms  can  be  usefully  di^layed  (Jan  et  al.  1989)  in  tabular  form  as 
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TABLE  1.  ORDER-OF-MAGNTTUDE  ESTIMATES  OF  AXIAL  AND  SECONDARY  FLOW 

COMPONENTS _ _ 


Axial  Flow  Equation 

Tonporal 

(unsteady) 

Convective 

Viscous 

Representative  terms 

dii 

Bu  Bu 

u  — ,  u  

dt 

Bx  By 

Order-of-magnitude 

(i>U 

1/2  V,U 

vU 

estimates 

Lg’  a 

Secondary  -  Flow  Equation 

Temporal 

Convective 

Centrifugal 

Viscous 

du 

d\)  dv 

Representative  terms 

Bt 

u  — ,  1)  — 

dx  dy 

R 

Order-of-magnitude 

utV, 

V,  u  vf 

Hi 

estimates 

s 

Lg  ’  a 

R 

The  transverse  secondary  flows  are  primarily  of  a  swirling  nature.  In  flie  pulmonary  airways,  it  is  fliey 
fliat  are  largely  reqx)nsible  for  the  mixing  of  chemical  toxicants  with  air.  Such  flows  are  particularly  sensitive 
to  bifurcation  angles  and  radius  of  curvature,  as  they  engender  centrifugal  forces  on  fluid  and  solid  particles 
entrained  in  such  flows.  Under  certain  circumstances  to  be  identified,  these  swirling  flows  may  be 
diaracterized  by  an  inherent  unsteadiness.  In  tiiis  case,  the  centrifugal  forces  caused  by  particles  following  a 
curved  trajectory  are  offset  by  temporal  accelerations,  such  tiiat 

Vs(o~U^/R 

where  Vs  is  the  secondary  velocity  scale,  U  an  axial  velocity  scale,  and  R  a  radius  of  curvature  of  tiie  particle 
trajectory  wifliin  tiie  airway.  Adopting  a  generation  length  L,  an  airway  radius  a  and  a  breatiiing  frequency 
0),  one  obtains  an  order-of-magnitude  estimate  of  tite  secondary  velocity 

Vs/U~U/(oR  =  Ua  xa/R 
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where 


U  =  L  co/2. 

From  the  ob^rvation  that  the  transverse  accelerations  VgU  are  much  larger  flian  cither  the  convective 
accelerations  or  fee  viscous  fcs'ces,  it  follows  feat,  for  a  typical  human  airway  bifurcation  wife  a/L  =  1/7,  and 
a/R  =  0. 1,  Womersley  munber  =  a*co/v  must  be  much  larger  than  one,  if  the  flow  is  to  be  dominated 
by  unsteadiness. 

Similar  order-of-magnitude  analyses  have  been  carried  out  by  Jan  et  al.  1989  for  fee  viscous  and 
convective-dominated  flow  regimes. 

Their  results  are  aptly  summarized  in  Figure  1  below: 


Figure  1.  Flow  regimes  based  on  order-of-magnitude 
analysis  mapped  onto  a  log  scale.  The  3  limiting 
regimes  are  I  (unsteady),  II  (viscous),  and  III 
(convective),  defined  in  terms  of  effects  that 
dominate  in  each  zone.  See  text  for  further 
explanation. 
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By  locating  the  flow  regime  applicable  to  given  animal  or  human  tests  within  the  context  of  the 
above  Figure  1,  it  is  straightforward  to  classify  those  flows  within  the  aforementioned  boundaries  as  being 
dominated  either  by  viscous,  convective  accelerative  or  unsteady  forces.  Extrapolation  of  such  test  results 
should  only  be  logically  carried  out  amongst  species  in  which  the  same  flow  regime  apphes. 
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